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An optical model of classical photons propagating through array of many beam splitters is devel- 
oped to give a physical analogy of Parrondo's game and Parrondo-Harmer- Abbott game. We showed 
both the two games are reasonable game without so-called game paradox and they are essentially 
the same. We designed the games with long-term memory on loop lattice and history-entangled 
game. The strong correlation between nearest two rounds of game can make the combination of two 
losing game win, lose or oscillate between win and loss. The periodic potential in Brownian ratchet 
is analogous to a long chain of beam splitters. The coupling between two neighboring potential 
wells is equivalent to two coupled beam splitters. This correspondence may help us to understand 
the anomalous motion of exceptional Brownian particles moving in the opposite direction to the 
majority. We designed the capital wave for a game by introducing correlations into independent 
capitals instead of sub-games. Playing entangled quantum states in many coupled classical games 
obey the same rules for manipulating quantum states in many body physics. 
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I. INTRODUCTION 

It is an old saying from Chinese philosophers 3,000 
B.C. that when a thing approaches to its extreme point 
step by step, it will finally transform into the opposite 
thing Scientific discoveries stem from different fields 
to expand this philosophical idea. In engineering me- 
chanics, combing two unstable system in the correct way 
may produce a stable system[2]. Parrondo devised a 
game to show the combination of two losing games may 
lead to an ultimate wining game 3]. This game paradox 
was used to explain the directional motion of molecule 
motors by combining two random Brownian motion ac- 
cording to Parrondo's game paradox}^. In quantum frus- 



trated spin system, quantum fluctuation or thermal fluc- 
tuation which is supposed to destroy ordering results in 
long-range magnetic ordering^, this is termed as order 
from disorder [6|. 

At first sight, Parrondo's game is counterintuitive. 
Suppose a player has capital 100 in beginning. He plays 
two losing game: game A and game B. If he only plays 
game A, the capital will drop. If he only plays game B, 
the capital drops too. If he plays game A many times, 
the capital drops, he switch to play game B many times, 
the capital continuous to drop, and then he switch to 
play game A again, and so on. One would intuitively ex- 
pect that the player will finally lose, if the player wins in- 
stead, the game may be called a paradox. In fact, if there 
is no any correlation between the nearest two rounds of 
game, no matter how the player switch between two los- 
ing games, the final game will lose, it is indeed a paradox 
if the game wins instead. But if there exist correlation 
between the nearest two rounds of losing games or any 
two rounds of game, the combination of two losing games 
could be a win, a loss or oscillation between win and loss. 
It depends how the correlation is planted among different 
games. 

Parrondo's game introduced the correlation between 
the nearest two rounds of game. Parrondo, Harmer and 
Abbott proposed a history-dependent game paradoxQ 
as an improvement of Parrondo's game[3(|. Later on, the 
history-dependent game paradox was extended to couple 
two history-dependent games 8]. The history-dependent 
game paradox[7j consists of two games: game A and game 
B. The game A is a biased coin which has probability p 
to win and probability 1 — p to lose. The players of game 
B are four coins whose strategy at each step depends on 
previous two steps in history [7]. Both Parrondo's game 
and this history dependent game are mathematical is- 
sue which study the relationship between the winning or 
losing probabilities of two coupled games. 

An indivisible classical photon scattered by a beam 
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splitter is a perfect analogy of coin. If one plays a coin 
100 times, it faces up 80 times and faces down 20 times, 
we say the winning probability is 0.8, and losing probabil- 
ity is 0.2. Each time the coin faces up, the point increase 
by +1, on the contrary, we has —1 point upon input cap- 
ital. The final gain of this game is (+80 — 20 = +60). 
Therefore this is winning game. We let a photon col- 
lide with the beam splitter 100 times, it passed the beam 
splitter 80 times and is reflected by the beam splitter 
for 20 times. The reflection coefficient is 0.2 and the 
transmission coefficient is 0.8. Photon is electromagnetic 
wave. As we know, there is a e"^ phase shift whenever 
a wave is reflected. Thus the reflected photon naturally 
carries a negative sign e^^ = — 1. When two waves meet 
each other with tt phase difference, the destructive inter- 
ference reduce the overlapping intensity to zero. In the 
language of photon, negative photon will annihilate with 
positive photon. If we let 100 uncorrelated photons col- 
lide with the beam splitter, and put the reflected photon 
and transmitted photon together, the number of survived 
photons is (+80 — 20 = +60). The classical photon is not 
only analogy of a coin, but also a solid physical imple- 
mentation. 

I will use photons passing through beam splitter ar- 
ray to give an almost exact implementation of both Par- 
rondo's game and the history-dependent game. The dif- 
ference is Parrondo-Harmer- Abbott game does not have 
the negative sign —1 for the lost coins in their mapping 
matrix. While I will keep the e" phase shift for every 
reflected photon. There are other ways to realize similar 
probability relation, such as special dices, quantum parti- 
cles passing two slits, an so on. A light beam propagating 
through a network of mirrors(or more accurately beam 
splitters) is a very familiar phenomena for most people. 
This optical system give us a clear picture to show how 
the probability flow transform from loss to win, and vice 
verse. 

Playing game will produce a time series recording the 
instantaneous value of output capital. We map every 
time series into an unique path in spatial degree of free- 
dom. One can combine many losing or winning games 
by drawing optical diagrams. A real optical system pro- 
vide a practical way to test those complex graphical de- 
sign. Parrondo's game only combines two games. What 
I demonstrated in this optical model is a diagrammatic 
method to combine many games (the number is at least 
larger than three). Every local game may be a loss or 
win. When many of these games are connected to form 
a complex network, one has to calculate to see the final 
output. The correlation between neighboring games play 
the key role in determining the final output. 

The coupling between two neighboring beam splitters 
provide a direct understanding to the directional motion 
of Brownian particles in periodic asymmetric potential. 
The combination of the two losing games in Parrondo's 
game paradox was used to explain Brownian ratchet 01 ■ 
In fact, what really matters is the strong coupling be- 
tween neighboring potential wells. Every local potential 



well behaves like single beam splitter. The backward 
probability is turned into forward probability through 
two continuous steps of reflection by two neighboring po- 
tential wells. The optical model suggest that Brownian 
particle in periodic asymmetric potential does not always 
move in one designed direction, they can oscillate back 
and forth, or move in the opposite direction to the de- 
signed one. 

The article is organized as follows: 

In section II, we take a classical photon traveling 
across the array of beam splitters as analogy of Parrondo- 
Harmer- Abbott game. A different transfer matrix from 
that of Parrondo-Harmer- Abbott game is derived by con- 
sidering the phase shift of reflection beams. When both 
sub-games lost, the final output of this optical game could 
be a win, a lose or oscillating between a win and a loss. 

In section III, Parrondo game is also expressed into 
my optical analogy. The modular operation in Parrondo 
game is equivalent to a special beam splitter. The opti- 
cal diagram has the same structure as that of Parrondo- 
Harmer- Abbott game. 

In section IV, we designed the fractal history tree of a 
generalized history dependent game with longer memory, 
and extend it to the loop lattice of many beam splitters. 
We calculated the numerical eigenvalues of the transfer 
matrix for typical cases. When all the sub-game lost, the 
combined game may has oscillating output, or win, or 
loss. 

In section V, a game with entangled history states is 
designed. For a symmetric state, two entangled winning 
game give a losing game. For antisymmetric states, the 
two games liken to lose in order to win the combined 
game. The final output of combined game oscillates be- 
tween loss and win if only one of the two can win. 

In section VI, it shows the optical analogy of Brownian 
ratchet and how to play a quantum wave using classical 
capitals. A general way of combining many quantum 
sub-games is proposed. It also discussed the relation be- 
tween quantum many body Hamiltonian and the com- 
bined game of many sub-games. 

The last section is summary. 



II. AN OPTICAL ANALOGY OF 
PARRONDO-HARMER-ABBOTT GAME 

The game A in Parrondo-Harmer- Abbott game para- 
dox is played by a biased coin. The coin has probability 
p to face up, and the probability of facing down is 1 — p. 
If it faces up, the capital increase by one, otherwise the 
capital decrease by one. Suppose we have capital C in 
the beginning, after many rounds of game A, the winning 
capital approaches to C x (+1) x p. On the losing side, 
the gain isCx(— l)x(l — p). The final gain of the game 
is G — C{2p — 1). Usually the capital is positive, if the 
capital is negative, we will get a negative value from the 
winning output, C x (+1) x p < 0, if C <{). While the 
losing side give a positive value, C x (—1) x (1 — p) > 
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0, if C < 0. In my optical analogy, the negative capital 
is represented by the photon which is reflected for odd 
times. 

I take an optical beam splitter as analogy of a biased 
coin. A common Beamsplitters is a half-silvered mirror. 
Only one part of the input beam can pass through the 
mirror, the other part is reflected. The thickness of the 
coating aluminum controls the ratio between transmit- 
ting light beam and reflecting light beam. If the light 
beam is weak enough so that there is only one photon 
arriving at the mirror within certain time interval. If the 
photon pass through the mirror, the transmitted photon 
number will increase by one, the transmission probability 
is p. If the photon is reflected, then the total number of 
photons in the reflected beam will increase by one, the 
reflection probability is 1 — p. A photon is also electro- 
magnetic wave, every time it is reflected, a phase shift 
of expfiTr] will be attached to the photon. Therefore the 
losing photon spontaneously flip a sign. If the time inter- 
val between two nearest sequent photons is long enough 
to erase their correlation, we get a random sequence of 
game A. In the end, the ratio between the total num- 
ber of transmitted photons and total number of reflected 
photons is — p). 

Game A can be equivalently implemented by a ran- 
dom sequence of uncorrelated photons. The photon of 
game A has no memory of history. We have two equiv- 
alent way to play game A with many photons. One way 
is preparing an identical mirror for each photon, and let 
the photon pass it only once(Fig. [T]). The other way is to 
let many photons pass one mirror(Fig. [Ij. In the former 
way, we sum up the number of transmitted photon and 
reflected photon after all the mirrors to get the transmis- 
sion intensity and reflection intensity. In the most ideal 
case, the value of the intensity of the two beams are the 
same as measured from the single mirror for many pho- 
tons. Uncorrelated many photon systems is essentially 
classical light source. 

To inject an ideal classical light source into a beam 
splitter is equivalent to playing game A independently 
many rounds. The number of the rounds of game is the 
number of photons in the light beam. The intensity of the 
light beam is the capital. The beam splitter is the biased 
coin. The photons transmitted through the mirror are 
winners, while the reflected photons are losers. Suppose 
the intensity of the input light(or in other words, the 
initial capital) is Iq. The total gain of those winning 
games is Iqp, while the loss is counted by the reflected 
photon, /o(l — p)(Fig. [IJ. To make game A as a losing 
game, the reflected photons must be more than those 
transmitted. The transmission probability is less than a 
half, (p=l/2-e, e > 0). 

The game B proposed by Parrondo, Harmer and 
Abbott [3] can be interpreted as a single Brownian parti- 
cle which determines its strategy for the third step based 
on the winning or losing states two steps earlier. This 
game can be simulated by computer in a direct way. 
However it is hard to find ideal correspondence in op- 
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FIG. 1: A beam splitter with more reflection than transmis- 
sion as an implementation of game A in Parrando game para- 
dox. The double-faced box with a dark side and a grey side is 
the beam splitter. Many uncorrelated photon passing through 
one beam splitter is equivalent to many single photon passing 
through many isolated beam splitter. 

tical analogy. 

We designed an array of beam splitter, and let single 
classical photon propagate across this array. Here the 
classical photon is real classical object like small metallic 
ball, even though I called it photon. Every time the clas- 
sical photon pass one beam splitter, it is either reflected 
or transmitted as one, it never split into two halves. We 
mark the path of transmitted photon by 1 for it is the 
winning path. The path of refiected photon is labeled 
by 0. First, we generate the two-step-earlier histori- 
cal states. The earliest beam is generated by a single 
beam splitter, it gives one reflection and one transmis- 
sion. Then we place one beam splitter in each way of the 
output of the earliest beam splitter to generate the second 
earlier beams. Each of the two new beam splitters will 
generate a new pair of transmitted beam and reflected 
beam. Finally we have four historical states recording the 
output beams two steps earlier, |0t_20t_i), |0t_2lt-i), 
|lt_20t_i), |lt_2lt~i). If the photon is positive in the 
beginning, the pair of |0t_20t_i) and |lt_2lt-i) has pos- 
itive value for the photon has been reflected even times. 
While the output of |0(_2lt-i) and |lt-20t-i) has neg- 
ative value since the photon is refiected odd times. In 
fact, The two-step history game give out a physical un- 
derstanding to the so-called game paradox. Suppose both 
the beam splitters at t — 2 and t—1 has a reflection prob- 
ability 0.8, the transmission probability is 0.2. The final 
gain of path |0t_20t_i) is a positive value, 0.64. The gain 
of |lt_2lt-i) is 0.04. The sum of these two positive path, 
0.68, is already larger than 0.5. Therefore, the three los- 
ing games of the beam splitters at i — 2 and i — 1 lead to 
a winning results. 

To implement the game B proposed by Parrondo, 
Harmer and Abbott [t'I , we place the designed beam split- 
ters in game B for each of the four history states. Each of 
the four states only interact with one special beam split- 
ter m\ at time t. Each m\ has a transmission probability 
of bi, and a refiection probability of (1 — hi). The games 
rules are summarized as following(Fig. [2]), 



m\\Qt-2Qt-i) 
m?|0t_2lt-i) 



m. 



^'lt-20i_l) 



m^ |lt_2lt-i) 



6i|lt)-He"(l-6i)|Ot), 
62|lt)-f e'"(I-62)|Ot), 
63|lt)+e"(I-63)|Ot), 
64|lt)+e"(l-54)|0t). 



(1) 



Every history states has an unique destination. For in- 
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stance, if the photon is reflected by the beam sphtters at 
time t — 2 and i — 1, it will meet beam splitter ml. The 
beam splitters at time t—2 and t—1 are only used for gen- 
erating the initial historical states of game B. After the 
first two- rounds of game B, the newly generated beams 
out of ml becomes the history for the next round. Game 
B can repeat itself on the four designed beam splitters, 
ml, i = 1,2,3,4. We denote the four historical states as 
a vector, 

= [|0t_20t-l),|0t_2lt-l),|lt-20t-l),|lt-2li-l)].(2) 

The transfer matrix map the historical states one step 
forward in the direction of time, 



l^t^l.t) = T+\lbt-2J-l)- 



(3) 



In this optical analogous system, we can directly read out 
the transfer matrix r+ from Fig. [51 



T+ 



h 








b2 





63 

e'^'h 
64 



(4) 



where bi is the reflection probability, bi = 1 — bi. T+ 
is similar to the matrix given by Parrondo, Harmer and 
Abbott [tII, but it is not mathematically equivalent to that 
matrix. Here we attached a phase shifter to every re- 
flection, can not be transformed into the matrix in 
Ref.Q through elementary matrix operation. It has dif- 
ferent eigenvalue and different eigenstates. We calculated 
the stationary states \ipsta) which is invariant under the 
operation of T+, {ipsta) = T~^\ipsta), it turns out to be 
IV'sta) = [0,0,0,0]. While the stationary solution of the 
matrix in Parrondo-Harmer- Abbott game has complex 
algebra [7]. 

The transfer operator can be physically carried out 
by properly connecting the output with input. The out- 
put of each beam splitter has fixed destination designed 
by the game rule. We showed the fiow chart for imple- 
menting the operator in the optical network(Fig. [5]). 
The output beam are fused into designed input beams. 
The photon is trapped in a cyclic optical network. To 
read out the information of the beam after many rounds, 
one may disconnect the output from the input. If one 
use one photon to play game B, one must play it many 
times to determine the distribution of probability out of 
each beam splitter. An equivalent way of playing game 
B is to input a collection of many uncorrelated photons 
one time, and measure the intensity of output beams. 

To get a game paradox in real physical system, we 
must choose the beam splitters with proper refiection 
and transmission to lose both game A and game B, and 
make sure the combination of game A and game B finally 
wins. The game A is a single beam splitter, if the trans- 
mission probability is less than reflection probability, i.e., 
p < 1/2, we lose game A. As for game B, Ref. [3] gave 
the value of four winning probability to lose the game B 
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FIG. 2; The optical flowchart for implementing the game B. 
The flow lines with arrow guides the flow of photons. The 
closed optical loop demonstrate how to realize the transfer 
matrix for repeating game B. 

in Parrondo-Harmer- Abbott game. However, the matrix 
in Parrondo-Harmer- Abbott game is not mathematically 
equivalent to the transfer matrix I used in the optical 
game B. So their results is not suitable for determining 
the transmission probability of the four beam splitters 
here. 

I will analyze the optical game B in detail and calculate 
the distribution of transmission probability for losing this 
game. First we make a fair history states through Bt-2 
and -Bt-i in Fig. [21 both of which has transmission prob- 
ability p = 1/2. The initial input of m^ is pt-2Pt-i = 
0.25. The input of is {l-pt-2)e" {l-pt-i)e''^ = 0.25. 
The input of m^ and mf are both (—0.25) due to the 
phase shifter e*'^. At ml, the photon has probability bi 
to pass and 1 — bi to reflect. The transmitted photon 
will go to mf_^i for it loses at -Bt_i and wins at ml. The 
reflected photon at ml will bent back to ml_^_i. At t + 1, 
the number of photon goes to beam splitter ml_^_i are 
7V,,(+i, (i = l,2,3,4). 



iVi,t+i = 0.25e*''(l-6i) -0.25e"(l- 

iV2,t+i = 0.256i - 0.2553, 

iV3,t+i = -0.25e'"(l-52) + 0.25e"(l 

Ni.t+i = -0.2562 + 0.2564. 



63), 



(5) 
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At time i-t-n, the number of photons to ml_^^ is calculated 
in the same way by operating the transfer matrix n times. 



.(6) 



[iVi_o, -/V2,0j A^3,0j A^4,o] is the initial number of photon at 
t = 0, here = [0.25, -0.25, -0.25, 0.25]. If the photon 
is reflected in the latest round, we call it a losing path, 
otherwise, it is winning path. At time t + n, the losing 
paths are \lt+n-i^t+n) and |Ot+„_iOt+„). The winning 
paths are \Qt+n-i'^t+n) and |lt+„_ilt+„). Therefore to- 
tal number of winning photons are N^m = A^i.t-Hn + 
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N4,t+n, the losing photons are Niose = A^2,t+n + -/Vs.t+n- 
If we want to win the game at t+n, then |-/Vi„i„| > |iV/ose|- 
We perform similarity transformation upon both sides 
of Eq. ([5]) to diagonalize the transfer matrix without 
modifying the states vector A''. Each eigenvalue is an 
algebra equation covered many pages. We denotes them 
as [Ai, A2, A3, A4]. The diagonal transfer matrix reads Eq. 
© greatly, 



(7) 



If the distribution of the transmission probability of four 
beam splitters satisfy the inequality. 
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KNi^o + A?A^4,o + A^7V2,o + A^TVa.o > 0, 



(8) 



the game wins at time t + n. This inequality equa- 
tion is defined within four dimensional cubic, {bi g 
[0,1], i = 1,2,3,4}. We divide this four dimensional cu- 
bic into 16 sub-cubic by dividing each bi into two halves, 
e (l/2,l],i = 1^2,3,4} and {b^"'" e [0,l/2),i = 
1,2,3,4}. One way to define a game paradox is the 
solution of inequity equation ([8]) exist in the sub-cubic 
|^;ose g [0, 1/2), i = 1, 2, 3, 4}. Notice that there are four 
unknown variables but only one inequality equation, the 
solution space is a three dimensional manifold. 

I showed four numerical eigenvalues of the transfer 
matrix in table HI /3k — ^2, &3, ^4} represents one 
distribution of the four transmission probability. We 
calculated four special cases: (3i = {bi — 0.2, 62 — 
0.2, 63 = 0.2, 64 = 0.2}, ^2 = {0.8,0.8,0.8,0.8}, 
^3 = {0.2, 0.2, 0.8, 0.8}, P4 = {0.8, 0.8, 0.2, 0.2}. The cor- 
responding eigenvalue of the transfer matrix is shown in 
table m Pi represents four losing beam splitters, the cor- 
responding eigenvalue Ai is negative, A2, A3 and A4 are 
almost zero. If we play game B odd times according to /3i, 
the game loses, while play it even times, the game wins. 
/3i is a solution of the inequality equation. /32 represents 
four winning beam splitters. It is also one solution of 
the inequality equation. jSs and /34 represents two losing 
beam splitter and two winning beam splitters. For 
Ai and A2 are imaginary numbers, if we play game even 
times, the contribution of Ai and A2 are negative, while 
A| and A4 are positive. But the initial input has nega- 
tive value for A2 and A3, thus the sum of the four paths 
are almost zero. The number 0.77 in the table is the 
approximation of 0.774597. There exist very subtle dif- 
ference between different eigenvalues. The above is just 
four examples of many solutions. I also checked numer- 
ically many other case, there are many solutions within 
the losing sub-cubic {6^°"^ e [0, 1/2), i = 1, 2, 3, 4} which 
can reach a winning result. 

In analogy of the combination of two games in 
Parrondo-Harmer- Abbott game|^, I designed the opti- 
cal network to combine the optical game A and optical 
game B here. It is almost the same network for imple- 
menting game B(Fig. [5]). The beam splitters used for 
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TABLE I: The four eigenvalues of the transfer matrix for 
different transmission probability distribution: /3i = {61 = 
0.2, ba = 0.2, 63 = 0.2, 64 = 0.2}, /Sa = {0.8,0.8,0.8,0.8}, 
Ps. = {0.2, 0.2, 0.8, 0.8}, Pi = {0.8, 0.8, 0.2, 0.2}. 




FIG. 3: The optical flowchart for playing the combined game 
in Parrondo-Harmer- Abbott history dependent game. The 
closed loop repeats the game sequence A^A-^B^B ■ ■ ■ . 



generating historical states for game B are now replaced 
by the beam splitter of game A. Eight beam splitters are 
added to play one more round of game B. When the single 
photon propagates through the array of beam splitters, 
it shows exactly repeating game A twice followed by two 
rounds of game B. 

First, the photon meet the beam splitter A at time 
t — 3(Fig. |3]), it has probability p to pass, and probability 
(1 — p) to be reflected. In either case, the output photon 
will be guided to meet another A-type beam splitter at 
t — 2. There will generate four possible states after the 
two A-type beam splitters at t — 2. Until now the pho- 
ton played game A twice to generate the initial historical 
states for playing game B. At time t — 1, we place the 
four beam splitters of game B in each way of the four 
historical paths. It will generate eight possible paths. 
Following the game rule between neighboring historical 
states in game B, we lead each of the eight possible paths 
into another eight assigned beam splitters at time t. So 
far we played both game A and game B twice. 

We input uncorrelated photons one after another into 
the optical network started with A-type beam splitter 
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at time t — 3, and counts the total number of photons 
running out of each B-type beam sphtter at time t. The 
total number of transmitted photon is what we win in 
the combined game. The total number of photons out of 
the reflection path is the loss. To repeat the cyclic game 
for n steps by single photon, 



-^t-3^t-2Bt-iBtAt+iAt+2 



Bf+n-lBt+n, (9) 



we collect the transmitting path out of the B-type beam 
splitters at time t into one ultimate path, and fuse all the 
reflected paths at time t into another ultimate path(Fig. 
[3]). The difference between the two ultimate paths is the 
net gain. Then we bent the net gain back to the A-type 
beam splitter at t — 3. The game, A ^ A ^ B ^ B, 
starts all over again (Fig. [3]). 

If we do not use single photon to play the game, an- 
other equivalent way is to input many uncorrelated pho- 
tons into beam splitter j4_3 at the same time. If the 
transmitted photons out of the eight beam splitters at 
t is more than those reflected, the game wins. In the 
language of classical optics, the intensity of transmitted 
light is stronger than the reflected light. In Fig. [3l we 
use a positive circle to denote the collection of transmit- 
ted light beams, and a negative square to represent the 
reflected beams. 

This optical game is not equivalent to the mathemat- 
ical game paradox proposed by Parrondo, Harmer and 
Abbott |7i|. But it can help us to understand those math- 
ematical issues. For physicist, there is no paradox in this 
optical network. The intensity of reflected light repre- 
sent exactly the losing probability. The transmitted light 
measures the winning probability. Even if there are more 
reflection than transmission at a local beam splitter, as 
part of the reflected beam can transform into transmit- 
ted beam through other beam splitters, it is not paradox 
if the transmitted beam out of the last round of beam 
splitters is stronger than the reflected beam. 

In fact, one may design much more complicate optical 
networks for more optical games. We first construct an 
ultimate optical network that produce stronger transmit- 
ted beam than reflected beam. Then we decompose the 
total network into small networks, and check if the small 
network generate stronger reflected light than transmit- 
ted light. Every beam splitter introduce a variable. For 
a large network with m beam splitters, the solution of 
inequality equation of the final gain is distributed in a 
m-dimensional cubic. As the dimension grows, the so- 
lution space also expands. One can find many different 
ways of decomposing the total network, and each sub- 
network give more reflection than transmission. 



III. THE OPTICAL MODEL IN ANALOGY OF 
PARRONDO'S GAME 

Parrondo's gameQ is the combination of two games. 
Game A is the same as that in the history-dependent 
game, it has probability p to win and probability 1 — p to 



lose. Game B flrst check wether the capital is divisible 
by a number M or not. If the capital can be divided by 
M , it play coin 1 which has a winning probability bi and 
a losing probability 1 — 6i, otherwise it play coin 2 which 
has probability 62 to win and probability 1 — 62 to lose 3. 

I will design an optical system to theoretically visual- 
ize Parrondo's game(Fig. 2]). First we represent every 
coin by a beam splitter. The capital is a collection of 
many classical photons. Game A is realizable by single 
beam splitter. Game B needs two different beam split- 
ters. The main difficulty of implementing Game B lies in 
doing the modular algorithm with photons. It is a very 
complex procedure to check if the number of photon is 
divisible by M or not. I just use a hexagonal box to rep- 
resent a fictional machine which can count the number 
of photons and do the modular computation(Fig. I 
call this machine a modular beam splitter for it behaves 
like an effective beam splitter. For example, if M = 5, an 
arbitrary number of input photons Ninput must be one 
element of the set {5fc, 5A: 1, 5fc 2, 5fc 3, 5fc 4}. We 
don't know wether Ninput can be divided by 5 or not. 
The probability for a divisible Ninput by 5 is 1/5. Ninput 
has a probability 4/5 of being unable to be divided by 
5. So the Ninput photons has probability 4/5 to meet 
beam splitter B2, and has probability 1/5 to meet Bi. 
If we input different number of photons into the mod- 
ular beam splitter for many times, and denote the sum 
of these input number as Nt. The total number of pho- 
tons collected by Bi is Nt/5. While 4iVt/5 of the total 
number goes to -62- In this sense, the modular beam 
splitter is equivalent to optical beam splitter. However 
the modular beam splitter does not add phase shifter on 
its output beams. An optical beam splitter add a on 
its reflected beam. 

We play this optical game B by bending the output of 
Bi and B2 back to the initial modular beam splitter, the 
input capital between the nearest two rounds of game B 
reads, 



Nn+l = Nn 



— -(2b2-l) + —(2bi-l] 



(10) 



Suppose the initial capital is iVoi then the capital after k 
rounds of game is. 



Nk = NQ 



— — ^-(262 - 1) + — (26i-l) 



(11) 



For the special value of the probability parameter (61, 62) 
and integer M in Ref. Q, (M = 3, 62 = 1/10, 
bi = 3/4), the capital at the fcth round of game B is 
Nk = 7Vo(-ll/30)'=. If k is odd number, Nk < 0, the 
game lost. The game wins for an even number of k. If 
both (61,62) are larger than 1/2, the net output capital 
Nk would be constantly positive. For an extreme case 
of 61 = 62 = 1, i-e., Bi and B2 has zero probability 
to lose, then Nk ~ iVo[-|-l]'^. For another extreme case, 
61 = 62 = 0, Nk = Nq[—1]'^ which is positive for k is even. 
This is because if we input a (—1) into a completely losing 
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FIG. 4: The optical flowchart for implementing Parrondo's 
game. The hexagon with M represents the modular beam 
splitter. The circles with numbers represents the label of 
beam splitter. Small grey window are the symbol of photon 
detector. 

game, the output would be +1. In fact, a losing game B 
has two branches of output. If the computer simulation 
assumes the input capital can only take positive value, 
the positive branches out of negative input will be elimi- 
nated. Only the negative branches is left. Then one will 
see the capital decreases monotonically. The simulation 
result of game B as showed in Parrondo's game is only 
the negative branch. 

Parrondo's gameQ is a winning game by periodically 
repeating game A twice followed by two rounds of game 
B. Fig. |4] shows the optical flowchart for playing the 
cyclic game oiB^B^A^A^B^B---. The 
output of game A is brought back to the input of game 
B, it forms a closed loop. The game may start from any 
point in the optical loop. Every transmitted beam out of 
one modular beam splitter carries a number 1/M. The 
weight of the reflected beam from modular beam splitter 
is {M—l)/M. The reflected beam out of Bi and B2 is at- 
tached by an additional phase shifter e''^ upon their cor- 
responding reflection probability, i.e., e*'^(l — Playing 
two rounds of game B from single modular beam splitter 
give out 16 different output path. The final output of 
each path is calculated by multiplying the weight num- 
bers on all the bonds along this path until it reaches the 
photon collector (the small grey window in Fig. 0]). The 
input for beam splitter A at the third round is the sum 
of all the 16 paths. Beam splitter A has reflection prob- 
ability p and transmission probability (1 — p)e*'^. After 
two rounds of game A, the output flows back to the ini- 
tial modular beam splitter of game B, then the combined 
game starts over again. 

We can play this game many times using single classi- 
cal photon, it will give the same results as we play it once 
using many uncorrelated photon. The game procedure is 
almost the same as that in the historical game. The main 
difference is the modular beam splitter. Mathematically 



t+l^I) 




FIG. 5: The fractal tree array of beam splitters the game with 
3-step memory. The and 1 on each bond are label of losing 
state and winning state, t — 3 is the origin point. As t goes 
to inflnity, the fractal tree expand larger. The diagram here 
is the cutoff picture until t + 1. 

there is no difference between dividing iV by M and di- 
viding 1 by M/N. All the analyze about reflection and 
transmission probability holds for single photon. If we 
do not divide the single photon into many pieces, it is 
either reflected or transmitted at modular beam splitter. 
The reflection probability is (M — 1)/M, and transmis- 
sion probability is 1/M. The modular beam splitter does 
not results in e" phase shift on its reflection. 

Generally speaking, Parrondo's game and Parrondo- 
Harmer- Abbott history dependent game are the same in 
my optical analogy, except the transfer matrix in my op- 
tical model is not mathematically equivalent to that in 
their mathematical game. From a physicist's point of 
view, both game are reasonable without paradox. 

IV. COMBINATION OF GAMES WITH 
LONG-TERM MEMORY 

We can extend the memory of the game B to n-steps of 
history, n > 3. Every strategy for the next step depends 
on the historical states n-steps earlier. So we need at 
least n beam splitters. Every beam coming out of one 
beam splitter at every step can be clearly marked by 
reflection beam or transmitted beam. We can use new 
beam splitter to divide every old beam into two beams — 
a transmitted beam and a reflected beam. As long as 
there is no overlap between any newly generated beams, 
we construct a fractal tree of beam splitters. I will show 
one extended game with 3-steps of historical memory in 
the following(Fig. [5]). 

Whenever a photon reach a beam splitter, it is either 
transmitted or reflected. There are eight possible states 
if the photon pass through three beam splitters one by 
one. Repeating single photon propagation many times is 
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equivalent to propagate many uncorrelated photons once. 
A classical Light beams is a collection of many uncorre- 
lated photons. We count the winners by measuring the 
number of passed photons. The losers are numerated by 
reflected photons. Every path from the beam splitter at 
i — 3 to a beam splitter at t — 1 ends up with an unique 
beam splitter marked by binary code sequence(Fig. [S]), 

It := |0t_30t-20(-i), 2t := |0(_30t-2lt-i), 

3t := |0t_3lt_20t_i), 4t := \Qt-3^t-2^t-i) , 

5t := \lt-3Ot-2Ot-1), Gt := |lt_30t_2lt-i), 

7t := |lt_3lt-20t_i), 8t := |lt_3lt-2lt-i). (12) 

We introduce a beam splitter operator, m^, to express 
the action of the beam splitter jt at t. Beam splitter jt 
let pj of the beam pass, and reflect 1 — pj of the beam, 

m^|st_3St_2St-i) =Pj|lt) +e"(l -pj)|Ot), 

j - st_3 • 2" + st_2 • 2^ + st_i • 2" + 1. (13) 

Here j corresponds to the number enclosed by a small 
circle in Fig. [5j Since every beam splitter corresponds to 
an unique 3-step historical states. For simplicity, we use 
the sequence of beam splitters at any time t to express 
the corresponding 3-step historical states before 

^t^[U, 2t, 3t, 4t, 5t, 6t, It, %tf ■ (14) 

To replay this 3-step history-dependent game, the beam 
coming out of every beam splitter at time t must be 
guided back into the correct beam splitter at i + 1 . The 
mapping operator from state "0* to state V't+i is 

>i P5 0' 

pi P5 

p2 P6 

+ ^ P2 P6 , . 

3 P3 p7 ' ^ ' 

P3 p7 

p4 P8 

_0 p4 p8. 

where pi is the reflection probability attached by a phase 
shifter, pi = e"(l — Pi). T+ is the transfer matrix be- 
tween the nearest neighboring time point. The states 
after n rounds of game is reached by operating on 
the initial states n times, 

|St+„_3St+„-.2St+„-l) (T+)"|St_3St_2St-l). (16) 

|st_3St_2St_i) represents the initial input generated by 
three historical beam splitters. Following the fractal-tree 
in Fig. [SJ every path must pass three beam splitters, 
{mt_3, mt_2, TOt_i}(here we use time index to label the 
same type of beam splitter for simplicity), to reach the 
beam splitter at time t. There are one mt_3 beam split- 
ter, two mt-2 and four mt-i beam splitters. The initial 
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TABLE II: The five special distribution of the eight beam 
splitter's transmission probability. 



intensity of the beam can be read out by direct inspection 
of Fig. El 

P3 = Pt-3Pt-2Pt-l^ Pi^Pt-^Pt-^Pt-l^ 
P5 = Pt-Jt-2Pt-l, Pe = Pt-3Pt-2Pt-l, 

P-7 ^ Pt-:,Pt-2Pt-i-, Ps ^ Pt-3Pt-2Pt-i- (17) 

To play a fair game, the initial capital of eight paths, i.e., 
the initial number of photons in eight paths, must be the 
same. Thus we assign the same transmission probability 
to the three historical beam splitter, i.e., {Pt~3 — Pt-2 ~ 
Pt-i = 1/2). 

If we investigate the states of game after finite num- 
ber of rounds, it is convenient to diagonalize the transfer 
matrix T:^ by performing similarity transformation upon 
Eq. p^ . Then (T^)"- is expressed by the eight eigenval- 
ues A^(^-l,2,3,'••• ,8), 

(r3+)" = Dtag[X^„ \-,X-, A^, A^, A^, A?, A^]. (18) 

I calculated some numerical value of these eigenvalues for 
five special cases pl| : (1) /3i, eight beam splitters all lose; 
(2) /32, eight beam splitters all win; (3) fS^, four beam 
splitters lose, the other four win; (4) /34, six beam split- 
ters lose, the other two win; (5) /35, two beam splitters 
lose, the other six win. The corresponding eigenvalue are 
shown in table Hill The zeroes in the box are approxima- 
tion of 10~^° ~ 10~^^. The eigenvalue for the five spe- 
cial cases all includes complex number, and neighboring 
complex eigenvalue are conjugate. A7 and Ag are approx- 
imately zero. We study the all-losing case /3i as an ex- 
ample. The imaginary part of are Ai and Ai are 88 times 
larger than their corresponding real part. We take their 
real part as zero. Then wc play the game 2 times with /3i, 
the total gain is calculated by G = X]i=i Pi^i — 0.05355, 
so the final game wins even if all the eight beam splitters 
loses. 

A more geometric way of visualizing this game is to 
draw a fractal tree of history. Every bond on the tree 
has a number representing reflection or transmission. We 
follow each path from the origin point to the end, and 
multiply all the numbers on the bond along this path, 
then we get the gain of this path. The sum of all different 
paths give the total gain. For example, there are many 
future paths starting with 4( at time t. These future 
paths are given by the transfer matrix . In Fig. [51 
we showed the fractal tree of the paths started from At 
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TABLE III: The eight corresponding eigenvalues with respect 
to the transmission probability distribution in table |ll] The 
zeroes in small boxes are approximation of 10"^" ~ lO"^'"*. 
Im indicate the imaginary part of the eigenvalue. 
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FIG. 6: The decimal labeling of all the paths starting from 
beam splitter 4 until time t + A. The small black disc indicates 
the reflection beam. The small white disc represent the trans- 
mission beam. The numbers within 1 ~ 8 after each black or 
white disc are the label of the beam splitter in game B. 

until t + A. Any path is marked by an unique binary code 
sequence. Every three nearest neighboring binary code 
indicates one historical states. For instance, the binary 
sequence 

Ot_3lt-2lt-lltOt+iOt+2lt+3 (19) 

is equivalent to a sequences of 3-step states, 

|0t_3,lt-2,lt-l) ^ |lt-2,lt-l,lt) ^ |lt_l,lt,Ot+i) 
^ |lt,0t+i,0t+2) ^ |0t+i,0t+2,lt+3>. (20) 

In mind of the unique decimal labeling of the beam split- 
ter after the eight historical states, Eq. (IT^ . we get 

the sequence of beam splitters with different transmis- 
sion and reflection. 



1 1+2 



H+3 



(21) 



Every sequence like the above corresponds to an unique 
path in the fractal-tree network(Fig. This path is 
labeled by 57 (Fig. ^ which is the decimal mapping of 
the sequence of binary code Eq. (ITOl) , 

• 2^ -f 1 • 2^ + 1 • 2"^ + 1 • 2^ + • 2^ + • 2^ + 1 • 2°. (22) 

For a given time i', the path ended up with 1(/ is win- 
ning path. The paths ended up with Of are losing paths. 



The gain(loss) of winning (losing) path is calculated by 
the scalar product of the probability at every beam split- 
ter along the path. We calculate the profit of the path 
sequence Eq. until i + 3 as an example. We denote 
the number of photons of at beam splitter 4t at time t as 
iVo. A'o is the ultimate result of the historical path be- 
fore t. Following the path lfOt+iOt_|_2lt+3, A^oP4 of the 
Nq photons will pass beam splitter At- So the gain of 
this path at time t is Gt = Nqp^. At t + 1, the photons 
will be reflected by the beam splitter ^t+i- The number 
of reflected photons at i + 1 is Gt+i = Gtps — NqPaPs, 
where ps = 1 — ps is the reflection probability of 84+1. 
Following the same calculation, the photons will be re- 
flected at 7t+2 and transmitted through 64+3. The final 
gain of the path 4t8t+i7t+25t+3 is 



G 



(57) 
t+3 



NoPiP8P7P5- 



(23) 



This is only the gain of path 57 until time t + 3. There 
are 2*" ~ 128 independent paths at t + 3, part of which 
is shown in Fig. [51 To find out wether the total game is 
winning or losing until time t + 3, we must calculate the 
profit of all the paths, and compare the ultimate number 
of transmitted photon and the number of reflected pho- 
tons. If the transmitted photon is more than the reflected 
ones, the game wins, otherwise the game loses. 

The total number of output paths grows following 2* 
as time goes on. If we investigate the behavior of the 
game in the long run, it is better to use differential equa- 
tions with continuous time. Every beam splitter at time 
t corresponds to an unique three-beam splitter sequence 
in the latest history. Thus if one knows the label of beam 
splitter, one can deduce the information about winning 
or losing in the latest three steps of game in history. 
For instance, if the label of four beam splitters at time t 
are {1, 3, 5, 7} , the game along the four paths at time 
t — 1 is losing. For the paths ended up with beam split- 
ters { 2, 4, 6, 8} at time t, the game at time t — 1 is 
winning. We summaries the eight beam splitters into a 
vector in Eq. (|14p. The transfer operator maps a 
vector ■04 to a new vector ■0s(i + 1) at t + 1. The new 
vector of eight beam splitters bear the same meaning as 
that at time t, while the time index increased one step 
further. Within a longer period, we take the time inter- 
val as continuous variable so that the difference between 
neighboring states in history can be denoted as a deriva- 
tive, "04-1-1 ~ 4't = dfipt- Then we map the equation of 
transfer operator into a differential equation. 



d_ 

di 



(24) 



where [T+ — /] play the role of Hamiltonian. / is the iden- 
tity matrix. A similarity transformation performed upon 
both sides of this differential equation can diagonalize the 
Hamiltonian operator by keeping the equation invariant. 
Properly choosing reflection and transmission probabil- 
ity of the eight beam splitter, one can get an invariant 
vector which keep the same value all the time. This static 



solution is the eigenfunction of equation T^tpt 
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We can design similar flowchart to implement the eight 
dimensional transfer operator as that for the two-step 
optical game B. The transfer operator redistribute the 
beams among different beam splitters following the dy- 
namic equation (|24p. The beams passing through some 
beam splitter may get stronger, while others may become 
weaker. The time dependent distribution is the solution 
of dynamic equation (j24p . Exactly diagonalizing the ma- 
trix [T+ — /], we derived the eigenvalues of the eight beam 
splitters, { Ej, j = 1, 2, 3, • • • , 8 } . Each eigenvalue is a 
heavy algebra equation of the eight unknown transmis- 
sion probability. The solution of the dynamic equation 
Eq. p4|) has the following form. 



j{t)^Pjexp[Ejt], J = 1,2,3,..,! 



(25) 



Pj is the initial intensity of light beams at time t = 0. 

We decompose the eigenvalues into real part and imag- 
inary part, Ej = Re.j + i Irrij. The real part determines 
wether the solution diverges or decay to zero as time 
goes to infinity. The imaginary part governs the oscil- 
lation of the probability of finding certain beam split- 
ter at time t. The fractal tree game in Fig. [5] shows 
{l{t+l),^t+l),b{t+l),7{t+l)) collect winning beams 
at time t-f 1, while {2(t -I- 1), 4(t -I- 1), 6(t -I- 1), 8(t 1)} 
absorb the losing beams. Thus the total gain of the game 
at time t is the sum of four winning paths — {1,3,5,7}, 

G(0 - E^:'^^P[^J*]' J = 1.3,5,7. (26) 



The total loss until time t is 



L{t) - ^Pfcexp[Sfct], fc = 2,4,6. 



(27) 



For a winning game, the gain is lager than the 
loss, \G{t)\> \L{t)\. For a losing game, |G(t)|< 
|i(t)|. The inequality equation only provide one con- 
straint on the eight unknown variables, i.e., G{t) :— 
G{t,pi,p2,P3,P4:,P5,P6,P7,Ps)- Thus there are many dif- 
ferent ways to choose proper transmission probability for 
each of the eight beam splitters. For example, we intro- 
duce some equations by hand. 



Re[Ej] > 0, Im[Ej] = 0, 
Re[Ek] < 0, Im[Ek]^0, 



j = 1,3,5,7, 

A: = 2, 4, 6, 8. (28) 



These equations makes the winning beams grow stronger, 
the losing beams gets weaker. The zero imaginary part 
eliminated the oscillation. 

The above procedure of designing optical network for 
game paradox has a straight extension to the n-step 
history-dependent game. However only classical game 
on network can be analyzed by the same strategy of an- 
alyzing the three-step history-dependent game paradox. 
By classical game, I mean every path has a definite state 
of wining or losing. If there exist closed path in net- 
work, two beams are allowed to inject into the same beam 
splitter from the opposite surface of the mirror, then the 



transmitted beams from one side may fuse into the re- 
flected beam from the other side. In that case, we have 
to draw a double line — one winning line and one losing 
line — to label the path. 

However, the beam splitter in the loop network does 
not have unique transmission probability which can be 
determined self-consistently by the same historical map- 
ping rule. For example, if we construct a closed square 
by placing a single beam splitter at t + 1 to meet the 
other three beam splitters 3t — >■ nit-i — 4t in Fig. \5\ it 
can not simultaneously satisfy the mapping rule of the 
vertical path and the horizontal path. The vertical path 
leads to 5t, the horizontal path leads to 7t- One way 
to keep the historical mapping rule is popping into the 
third spatial dimension: place 5t on top of 7t at the same 
lattice site of two dimensions. 

In fact, we can design more general game paradox by 
abandoning the historical mapping rule. First, we con- 
struct a network of beam splitters, and define its output 
channel and input channel. Second, the winning capital 
is the sum of all the transmitted beams. The losing is the 
sum of all reflected channels. Finally, we find the proper 
transmission probability of each beam splitter so that the 
total winning capital is larger than the losing channel. A 
special exemplar design is shown in Fig. [7] (a) which in- 
cludes only four beam splitters. Each beam splitter j has 
the probability bj to win and probability 1 — bj to lose, 
j = 0, 1, 2, 3. We start from game 0, every time the games 
loses, we switch to next game, 0— >-l— >-2— >-3— )-0---. 
The winning beam are collected by the output rhombus. 
Even if each beam splitter has very small transmission 
probability, (bj = ej <^1, j = 0, 1, 2, 3, ) after a enough 
long time, all the photons will be collected by the output. 
This is winning game in the end. 
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FIG. 7: (a) A loops game with four beam splitters. All reflec- 
tion beams are transfer into transmission beam, (b) The loop 
game on square lattice. Each bond is the overlap of many 
transmission and reflection. 

The history tree for playing many rounds of the square 
game with four beam splitters in Fig. [7] (a) also have a 
fractal structure. As long as we can distinguish the re- 
flected beam with the transmitted at each round of game, 
every path in history has an unique label. Whenever the 
beam flows into a beam splitter at some lattice site at 
certain time, the beam splitter must make another copy 
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FIG. 8: The history-entangled game with two beam splitter 
in game B. The grey bar represents beam splitters. Three A- 
type beam splitters are used for generating initial historical 
states. 

of itself and split into two different paths to separate the 
overlapped reflection beam and transmission beam. The 
history tree of the game becomes a huge fractal network, 
especially when we extend this square game into a square 
lattice of many beam splitters. 

In mind of the fact we are dealing with the beam of 
probability flow instead of physical beams so far, we can 
design a game of real light beam on square lattice of 
many beam splitters like Fig. [7](b). In that case, we can 
not distinguish reflection beam from transmission beam 
passing through each bond. Firstly, one has to define the 
output beam splitter and input beam splitter along the 
boundary. The we inject N photons into the input beam 
splitter, and collect the output photon. If the output 
number of photons is larger than N/2, the game is win- 
ning, otherwise, the game is lost. The lattice structure 
enclosed by the boundary may be very complex. To cre- 
ate a physical game paradox, we cut the whole network 
into two or more sub-networks in a proper way so that 
there are less transmitted photons than reflected photons 
through each sub-network. But as a whole network, the 
transmitted photon is more than those refiected ones. 

V. HISTORY-ENTANGLED GAME 

If the game ask the player to choose its strategy for the 
next step according to historical record, but the player 
lost some historical information. In this case, the player 
will select all possible historical states to determine its 
strategy. The final output the game not only depends 
on historical states, but also depends how the player or- 
ganize those states. I call this kind of game a history- 
entangled game. 

I take a Brownian particle with incomplete memory 
as player. The particle first travel across three identi- 
cal A-type beam splitter to generate four initial history 
states(Fig. |5]). The transmission probability of A-type 
beam splitter is p^. Its reflection probability is (1— p^). 
If we have capital N in the beginning, then the output 
capital of the four beams after A-beam splitters are(Fig. 

Coo = iV(l-p^)^ Cn^Npl, 

Coi = Cio^N{pI~p^). (29) 

To make a fair game without dependence on initial con- 



dition, we choose exactly p^ =1/2 for symmetric states, 
then Coo = Cii = 0.25, Cqi = Cqi = -0.25. 

When the particle meet beam splitter Bi and B2 at 
the third step, it only remembers wether the past two 
rounds of game has the same result or different result. 
For instance, if it remembers that it won only once during 
the latest two rounds of game, but it does not remember 
in which round it wins. If it wins or loses two rounds, 
the only memory kept in its mind is the results of the 
two rounds are the same. The strategy of the particle is 
to put the output of \Ot-2Ot-1) and |lt_2lt-i) together, 
and submit them to Bi beam splitter. The output of 
|0t_2lt-i) and |lt_20t_i) are combined into B2 beam 
splitter(Fig. H]). 

There are two different ways to mix two historical 
states. One way is adding them up so that switching 
the two beams arrives at the same state. We called it 
symmetric state, 

1(^1,*) = |0t_20t_l) + |lt-2lt-l), 

= |0t_2lt-i) + |lt-20t_i). (30) 

The other way is subtracting one beam from another 
beam. If we switch the position of the two states, it 
will generate a (—1) sign on the original state. We call it 
antisymmetric state, 

\lpl,t) = |0t_20t_l) - |lt-2lt-l), 

\lp2.t) = |0t_2lt-l) - |lt-20t-l). (31) 

Both the symmetric state and antisymmetric state are 
similar to the well known Bell state in quantum mechan- 
ics. However I just use Dirac bracket to denote state for 
convenience, there is no quantum state here. The sym- 
metric states is comprehensible from the point view of 
single Brownian particle. It is also convenient for com- 
puter simulation. Since both |0t-20t-i) and |lt_2lt-i) 
will meet Bi, it is natural to add up their probability. 
The anti-symmetric states can be simulated by computer, 
but it is hard to explain physically by probability theory 
of conventional Brownian particle. We keep the antisym- 
metric state as a choice of particle's strategy. 

We first study the long term behavior of the symmetric 
state in by repeating game B. Bi take \(f>i,t), let bi of the 
input pass and refiect (1 — 61). \4>2,t) will meet beam 
splitter B2- There also generate a pair of transmission 
and reflection at B2- When the four output of Bi and 
B2 at t go to the next round game at another pair of Bi 
and B2 at t -I- 1, the particle began to take the output 
at t as historical states to reorganize the latest 2-step 
history following the same game rule. It view |0i,t) as a 
losing state |0t), and |02,t) is winning state |lt)(Fig. [5]). 
Repeating the pair of Bi and B2 will generate a sequence 
of the pair of states, 

|0t) = [|</'i,t),|</'2,t)]^. (32) 
The transfer matrix mapping every state vector \(j)t) one 
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step forward into |(/>t+i) is 



e*'^(l - 61) 
bi e'' 



62 

(1-52) 



(33) 



If game B is played so many times that the hfetime of 
the whole game is almost infinite comparing with single 
round of game, we take state function as a continuous dis- 
tribution, and write the difference equation as the deriva- 
tive of l^t) with respect to time, \(t)t+i) — \(t>t) = dtlcpt}- 
Then the discrete equation, \4>t+i) = T^\4't), is equiv- 
alently mapped into the differential equation, dt\(f>) — 
[T^ — As the probability element in the transfer 

matrix have no dependence on time, we can diagonalize 
the composite matrix [T^ — I] directly to determine the 
long term behavior of the states. The two eigenvalues of 
the composite matrix are Ei — — 2, i?2 — —2 + 61 + 62, so 
the solution of this differential equation is. 

=Poexp[-2t], 102) -Foexp[(-2 + 61 + 62)^1.(34) 

Pq is the initial capital at time t = 0. The total gain at 
time t is 



G{t) = Po(26i - l)e-^* - Po(262 - l)e 



^-2+61+62)* 



,(35) 



Here the absolute value of initial input capital is Pq — 

0. 25. The initial input is determined by Coo = C'li = 0.25 
and Coi = Cqi = —0.25. The total gain G{t) at time 
t = 6 is plotted in Fig. ^ (a) . As time goes to infinity, the 
total gain approaches to zero, limt^oo G{t) — 0. When 
both game Bi and B2 wins, the combined game loses. 
The gain function in Fig. IH] (a) drops to negative in the 
region 62 > 1/2, 61 > 1/2. When 62 < 1/2, the total 
gain is positive and increasing as 61 increases from to 

1. Thus if Bi and B2 loses, the game can still win. If 
beam splitter Bi wins at 61 = 1, while 62 increase from 
to 0.32, the gain is positive and reaches its maximal 
point at 62 — 0.32. As 62 continuous to increase, the 
gain drops to zero and finally becomes negative. 



equivalent to playing the two rounds of game only once 
using two separate Brownian particle. This indeed is true 
for the symmetric states. But to play the antisymmet- 
ric states, there must exist a correlation between the two 
Brownian particles. Suppose two Brownian particle meet 
beam splitters At-2 and At-i, one wins twice, the other 
lose twice. When they get into Bi, either the double 
winning particle or the double losing one must flip a sign 
to its capital. One particle must check the sign of the 
other one to determine its own sign. If play the antisym- 
metric state using single Brownian particle, there exist 
a historical correlation between the two-winning-step at 
one time and the two-losing-step at another time. If the 
single Brownian particle wins twice, when it meets Bi, 
it can not decide its own sign. So it has to search the 
history record for the nearest two-losing-step to check its 
sign. The computer can flip a sign of a number in an easy 
way. In physics, it is a complex design to add a exp[i7r] 
phase shift on a real photon, besides this phase shift can 
only be added after it transmitted through two nearest 
neighboring beam splitters. One may use optical medium 
to delay the photon. However, if we play many rounds of 
game, the number of optical medium will exponentially 
increase to infinity. 

The long term behavior of the antisymmetric state is 
significantly different from symmetric states. The game 
rule of beam splitter Bi and B2 for antisymmetric states 
is the same as that for symmetric states, 



B2\ 



6i|lt)+e"(l-6i)|Ot) 
62|lt)+e"(l- 62)10*) 



(36) 



The equation of motion for the antisymmetric state vec- 
tor, IV^) := [ 1^1), 1^2) r, is = [T+ - 1^. The 

transfer matrix between two nearest neighboring histor- 
ical states is 



T+ - 
J-i, - 



e"(l-6i) -62 

61 -6^(1-62) 



(37) 




[T^ — I] may be view as an effective Hamiltonian which 
has two eigenvalues, 

E- = i(-2-62+6i- VA), 

E+ = i(-2-62+6i + \/A). (38) 



(a) (b) 

FIG. 9; (a) The final gain of symmetric states at t—6. Here 
we chose Pq = 1 in the computation for convenience, (b) The 
final gain of antisymmetric states at time t=10. 

The output of antisymmetric states is much more com- 
plex than the output of the symmetric states. If two 
Brownian particle has no correlation, playing the two 
rounds of game twice using one Brownian particle is 



where A = (2 + 62 — 61)^ — 862 is the gap between two 
levels. Both the two eigenvalues are negative(Fig. [TUl) . 
So the two antisymmetric states will converge to zero 
when time goes to infinity. This a draw game in the 
long run. But its output in finite time domain could be 
winning, losing, or oscillating. We analyze the final gain. 

Git) = (26i - l)|^i) - (262 - 1)1^2), (39) 

from the time-dependent solution of the two antisymmet- 
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ric states, 

= PoCxp[i(-2-fe2 + &i-\/A)t], 

|i^2> = PoCxp[i(-2-fe2 + fei + \/A)t]. (40) 

The gain funetion G{t) has three different branches ac- 
cording to the gap A. 

(1) If 62 = 2 - 2^2 V&i + 61, then A = 0, two states 
become degenerated. In the degenerated states, if bi > 
&2, the game wins, otherwise the game loses. Notice that 
the results only depends on relative ratio between bi and 
62, even both bi and 62 loses, i.e., 61 < 1/2 and 62 < 1/2, 
the game can still be winning as long as bi > 62- On the 
other hand, if both 61 > 1/2 and 62 > 1/2, the game can 
also loses if 61 < 62. There also exist others case, such as 
one of the two game loses, the total game loses or wins. 

(2) If &2 < 2 - 2\/2^/Fl" + 5i, the gap is positive A > 0. 
The game 62 has to lose to make a winning combined 
game. In the winning region, when bi increases from 
to 1, the gain decreases(Fig. [S]). Thus 61 also inclined to 
lose to win the combined game. 

(3) If 62 > 2-2%/2V^+&i, the gap is negative, A < 0. 
In this region, the two game can not lose simultaneously. 
When bi loses, 62 is always larger than 1/2. If 62 lose, 
bi slide into the winning region. There is only one point 
of draw game, 61 = &2 = 1/2. However, this region 
allows the two game to win simultaneously, while the 
final gain of the combined game does not always win. 
G{t) oscillates between winning and losing(Fig. [TUl fb)). 
Although the gain is not real number but a imaginary 
number, in some sense, we still call this a game paradox 
since combining two winning game give a losing result. 




(a) " (b) 

FIG. 10: (a) The two eigenvalues of the transfer matrix for 
antisymmetric states, iJ_ and E+. They are both negative, 
(b) The final gain of antisymmetric states at t=17 within the 
oscillating region. 

The different behavior of the final gain between sym- 
metric states and antisymmetric states implies antisym- 
metric states is much more efficient to win the global 
game by losing every local games. For a history depen- 
dent game with long term memory, similar antisymmet- 
ric state can be constructed following the same strategy. 
The dimension of probability variables will grows higher. 
There is much freedom for one to design various different 
game paradox. 




FIG. 11: (a) The periodic potential is turned on and off. (b) 
The periodic potential is constantly presented. While the 
Brownian particles as a whole is in and out of the potential 
periodically or randomly. 

VI. OPTICAL GAME AS ANALOGY OF 
BROWNIAN RATCHETS AND MANY BODY 
PHYSICS 

A classical Brownian particle trapped in a periodic po- 
tential shows random motion, no matter the local po- 
tential well is symmetric or asymmetric. If the periodic 
potential with asymmetric local potential is turned on 
and off periodically, (Fig. HT] (a)), the Brownian particle 
would move in one direction Q. An equivalent way to see 
this phenomena is to fix the periodic potential but let 
all the Brownian particles together approach to the po- 
tential and leave there quickly (Fig. [11] (b)). This point 
of view describes the same thing by placing the static 
observer of reference coordinates on the potential well. 
Thus Brownian particles move to the right too. 

Every asymmetric potential well is analogous to a 
beam splitter. When the particle hit the potential well, it 
tends to stay in a lower potential point. If it collides with 
bond b in Fig. [ll](b), moving to the right is the dominant 
motion. When it hits bond a, it moves to the left(Fig. [TT] 
(b)). Switching the potential on and off in a periodic way 
in fact is equivalent to generate a standing wave of Brow- 
nian particles upon a static potential. In one dimension, 
Brownian particle has spontaneous random motion along 
the length of the periodic potential. If we ignore the colli- 
sion between different Brownian particles, one Brownian 
particle represents one propagating wave. As a rough 
phenomenological description, if a Brownian particle can 
cover an average distance \Brownian during one period 
of oscillation, AT, we define \Brownian as random wave- 
length. In one dimension, XBrownian is along the length 
of potential. The one domensional periodic potential is 
composed of many identical local potential wells which 
has unit length L(Fig. [TTjfb)'). 

If the random wave length XBrownian is smaller than L, 
the propagating wave only collide with one beam splitter, 
or one potential well(Fig. I12p . The transmission prob- 
ability is p = b/ {a + b), while the refiection probability 
is 1 — p = a/{a + b). Here the definition of transmission 
and reflection depends on the direction of input wave, 
one may switch a and b for different cases. If the ran- 
dom wave length XBrownian is larger than L, the prop- 
agating wave will cover two potential wells at the same 
time. In that case, the wave will pass two beam split- 
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FIG. 12: (a) The Brownian wave length within one period is 
smaller than the length of one local potential well, (b) The 
Brownian wave length is larger than the length of an unit of 
potential well. 

ters within one period. We must sum up all possible 
path within two-steps, it has been shown in my optical 
model of Parrondo-Harmer- Abbott game. But here the 
case is simpler for all the beam splitters are identical. 
If the random wave length XBrownian is so large that it 
covers many potential wells at the same time, we must 
take into account of three or more beam splitters at one 
time. This case meet the history dependent game with 
long term memory. 

The two or more steps of collision between neighbor- 
ing potential wells can be mapped into a limited num- 
ber of elementary optical diagrams. I showed the opti- 
cal diagram of two steps of collision in Fig. [131 The 
transmission is in the same direction as the incoming 
velocity. While the reflection is in the opposite direc- 
tion. Every reflection bond is labeled by a number, 
— {1 — p) = 6*^^(1 — p). If Brownian particle pass two po- 
tential wells continuously, we label it as state The 
final output particle moves in the same direction as the 
input velocity. The state |10) denotes that the Brownian 
particle first pass one potential well and then is reflected 
by the next potential well. The final output velocity of 
state 1 10) is in the opposite direction of the input ve- 
locity. So does the state |01). For the state |00), the 
Brownian particle is firstly reflected into the opposite di- 
rection at the first potential well, then it is refiected back 
to the original direction by the second potential well. The 
final velocity is in the same direction as input velocity. 
With this correspondence between optical diagram and 
collision within two neighboring potential wells, one can 
calculate the final output velocity of Brownian particles 
in flashing ratchet. Adding up the probability of |00) 
and 1 11) gives the positive velocity, 

F+ = P(|00))+P(|11)). (41) 

The negative velocity is determined by |00) and |11), 

F-=F(|01))+P(|10)). (42) 

If > \V~\, the particle moves in the same direction 
as its original velocity. On the contrary case, the particle 
move in the opposite direction. If \ V^\ = \V~\, the parti- 
cle is trapped in a local potential well and moves around 
randomly. If the length of unit potential well is much 
smaller than the phenomenological random wavelength, 
we must include the three-steps of collision or even more 
steps. The strategy is the same as I showed above. 




FIG. 13: The four elementary optical diagrams for a particle 
propagating through two beam splitters. The right hand side 
is the corresponding behavior of a Brownian particle colliding 
with two neighboring potential wells. 

If Brownian particles do not collide with each other, 
the physics of single particle is an ideal meanfield approx- 
imation of many particle system. In fact, the random 
wave length XBrownian is only well defined for a crowd 
of many Brownian particle. When many Brownian par- 
ticles are trapped in a periodic potential, the collision 
between different particles are modulated to form quasi- 
regular density wave. When this density wave propa- 
gates across the periodic potential, it is either reflected 
or transmitted. Modeling many colliding Brownian par- 
ticles as density wave is convenient for describing and cal- 
culating Brownian rachet in two dimensional periodic po- 
tential(Fig. [14] (a)). In two dimensional periodic poten- 
tial, a Brownian particle can circumvent potential wells 
instead of being blocked there. If the Brownian particle 
can randomly fluctuate around a circle, such as sperm or 
bacteria, the angular velocity is equivalent to an effec- 
tive magnetic fleld. In quantum Hall effect of electrons 
in magnetic field, a series of plateaus appear in the Hall 
resistance measurement. This give us a hint that the two 
dimensional flashing rachet for circling Brownian parti- 
cle maybe has signiflcant different phenomena from one 
dimension. 




(a) (b) 

FIG. 14: (a) The asymmetric potential well for Brownian 
rachet periodically distributed in two dimensions, (b) Ele- 
mentary flipper and its combination. Each box represents a 
flipper which corresponds to one losing game. A losing game 
is operated by odd number of elementary flipper. A winning 
game can be performed by even number of flipper. 

When the local density of Brownian particle oscillates 
periodically under the modulation of external potential, 
it may be called as density wave. It is hard to control 
exactly the local density, but in a game, a player can 
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determine the amount of his capital for each round of 
game. If we play the capital in a proper way, it will gen- 
erate a good correspondence between Brownian density 
wave and capital wave of game. 

First, we decompose the input capital as the inte- 
gration of a pair of conjugated complex function, N — 
J ip*ip = ipt is a wave function whose amplitude is 

the square root of the total capital, 

4' = Mt) + ^Mt) = V^exp[zl7(i)], (43) 

The angular velocity VL{t) = arctan[i/)f,(t)/?/'o(i)] is a com- 
plicate function of time. The capital now has inter- 
nal degree of freedom. Suppose we input capital 1000 
points within time interval T. But they do not join 
the game one after another in constant speed or once 
and for all. We divide those time interval into n pieces, 
T = [0,Ti) U [ri,T2) U [r2,T3) • • • U [T„_i,T„]. The in- 
put capital in each step is Nt^ which is oscillating within 
T, and = 1000. For example, we input an cosin 

wave, N = cos[fi(<)], t e T. Then Nt, = cos[Q.(t)]. 
The discrete capitals within T now are not unrelated 
numbers, they organized into waves. This give us the 
real part of the capital wave i?e[?/'] = ^/N cos[51(i)]. The 
imaginary wave of the capital is the delayed or advanced 
real wave by a phase of 7r/2, Imlip] ~ \/N cos[17(f)±7r/2]. 
The total input capital is 

N ^ Re[il:f + Im[il}f . (44) 

When a beam splitter received the capital Re[il)Y and 
Im[4'Y , it take them as independent beams and redis- 
tribute them separately. Both the real wave and delayed 
real wave(or imaginary wave) generate a transmission 
wave and reflection wave, we sum the transmission up 
and denote it as ■01. The sum of refection is tpo. The to- 
tal capital wave if the sum of the two waves, ip — tpi+ipo- 
We use Dirac bracket to express the winning wave as 
|1) Vi, the losing wave is |0) -00, = |1) + |0). 
Then the absolute value of total capital reads 

iV=(l|l) + (0|0). (45) 

We have assumed there is no overlap between reflection 
and transmission, thus (0|1) = (1|0) = 0. (1|1) = pii is 
transmission probability determined by the beam split- 
ter. (0|0) = Poo is the reflection probability. 

If the input angle between capital wave and beam split- 
ter is so small that reflection overlaps transmission, the 
transformation between reflection and transmission is not 
zero, i.e., (0|1) 7^ (1|0) 7^ 0. In that case, the total 
transmission is A'^i = + (1|0). The total reflection 

is iVo = (0|0) + (0|1). The reflected capital must carry a 
phase factor e^^ , 

7Vo = e"[(0|0) + (0|l)]. (46) 

A winning game requires more transmitted capitals than 
the reflected ones. The final gain of the game, G{t), is 

G{t) =Ni+No>0. (47) 



In the meantime, the sum of the two absolute values are 
the total number, N = |A^i| -I- |iVo|. We can define nor- 
malized transmission p — |A^i|/A^ and normalized reflec- 
tion 1 — p — |A''2|/A^ to get rid of the dependence on the 
absolute value of input. If the gain G{t) < 0, the game 
loses. 

Any game has an input and output. The output is 
either win or loss. If we input a positive capital, a pos- 
itive output corresponds to the win, a negative output 
corresponds to the lose. If we add a negative sign on 
the input, the output will flip a sign correspondingly. If 
the game either transform the positive input capital into 
a negative output or transform the negative input into 
positive output, this is always a losing game, we call the 
game a flipper. A flipper check the sign of input and sent 
out an output with opposite sign. The simplest flipper is 
a mirror with total reflection. 

Single game consisting of odd number of flipper is los- 
ing game, while the combination of even number of flip- 
per is a winning game. For example, we take a sequence 
of flippers, {A,B,C,D,- ■ ■} (Fig. [11(b)). The flipper 
A maps the input capital Ci„ to Cout under the con- 
straint Cin'xCout = —1, SO do the flippers {B, C,D,- ■ ■}. 
Different flipper produce different ratio between output 
capital and input capital, we denote this ratio Cout /Cm 
as ampliflcation, = C^ut/Cin^ ^^re {a — A,B,C,- ■ ■} 
is the label of the flipper. The ampliflcation of a flip- 
per is always negative. Fa < 0. If we play flipper A, it 
always lose, any positive capital turns out to be nega- 
tive. But flipper A also transform negative capital into 
positive capital. If we play two flipper A, the output of 
the first A is the input of the second A, then the posi- 
tive input capital at the first A will give out a positive 
output at the second A. The amplification of the double 
A game is Faa = FaFa > 0, Fa < 0. If we play a 
sequence of flipper {A,B,C}, the final amplification is 
Fabc = FaFbFc < 0. For the most general case, the 
product of a sequence of flipper, {ai, i = 1, 2, 3, • • • }, the 
final amplification is 

n 

F = Y[Fa^. (48) 

i=l 

If n = 2k + l, F < 0, it is losing game. If n = 2fc, F > 0, 
it is winning game. 

The input and output of any flipper must be classi- 
cal variables. The player is actually a classical observer. 
The player make decisions following classic logic. A flip- 
per may have complex internal game process dealing with 
complex quantum states. When we try to check the out- 
put of a flipper to see it wins or loses, we must mea- 
sure the weight of those states, all the results will be 
summarized into a classical number. Therefore the inter- 
junction between independent flippers is always classical 
operation, no quantum states could survival. Quantum 
states may be used to design the internal gaming process 
within a flipper. If two independent flipper are perfectly 
organized into one combined flipper, then the old inter- 
face between single flippers disappeared. The quantum 
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FIG. 15: The dark star in the center is where the game begins, 
input capital is injected there. The box with label 1 represents 
the direct product of odd number of elementary flippers. The 
box labeled by 2 represents the direct combination of even 
number of flippers. The cups attached to the wall represent 
quantum state collector. 



states can flow through those old classical channels. If we 
cut the united flipper into two parts without measuring 
the cutting edge, the quantum states are still waiting on 
the interface. But we can not determine wether each part 
wins or lose. They will collapse into classical numbers if 
we perform measurement. 

Strictly speaking, there is no quantum game paradox 
as long as we want a definite output state of either win- 
ning or losing for each sub-flipper. Quantum states can 
be used to design the internal structure of single flipper. 
For example, we choose electron with spin as player. The 
winning state is more spin-up than spin down, the losing 
state is the opposite case. The input electrons may be 
all spin-ups. Then we design various complex scattering 
process within the flipper, such as two slits, many slits, 
magnetic fleld, electric field, other fixed nuclear spin, and 
so on. On the output side, we measure the state of all 
output electron. If the total number of spin-ups is more 
than those spin down, the game wins, this is not success- 
ful flipper. This design must be abandoned and start a 
new design until we get more spin-down than spin-up. 
Once we get a flipper with quantum internal structure, 
and make a good interface between the output of one flip- 
per and the input another flipper, one can play the game 
paradox using many copies of the same flipper. From a 
global point of view, this is a classical game. If we go 
inside every flipper, it is a quantum processor. 

A single flipper does not results in so-called game para- 
dox. It is only when we combine many flippers in a com- 
plex way, the player can not read out the flnal output 
at first sight, the total game will produce seemly strange 
results. To find out the output of flipper network has a 
lot of things in common with many body physics. Espe- 
cially when the flipper are quantum processors, we are es- 
sentially manipulating various entangled quantum states 



from local flipper to reach a global winning output. I 
showed an example of combinate game on 5 x 5 square 
lattice(Fig. [T5|) . This game is the combination of 25 
sub-games. There are 14 losing games, and 11 winning 
games. The winning game are performed by the winning 
processor which is produced by even number of elemen- 
tary flippers. They are represented by the boxes labeled 
by 2. The processor produced by odd number of elemen- 
tary flipper are in charge of the losing games, they are 
labeled by i(Fig. [T5|) . The game starts from the dark 
star in the center. The dark star is a losing processor, it 
gives out four outputs. At the next step, the four out- 
puts flow into a ring of eight processors. Three of them 
wins(the disc with label 2 in Fig. [TS]) The other five would 
lose(the disc with label 1 in Fig. [15]). The 12 outputs 
of this ring provide the input for the last round of game. 
There are 16 processors in the last round, 8 lose and 8 
win. The output result at the first round of game can be 
read out by performing measurement on each of the 12 
output. The output states from winning processors are 
projected to final winning state (1|, this give the positive 
gain. While the output of those losing processors counts 
the negative gain. The total gain is the sum of all the 
outputs of this round, 
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the indexes i,i runs over all the eight processors repre- 
sented by disc in Fig. [151 If G{t) > 0, the game wins, 
otherwise the games loses. For a game combined n x n 
subgames, the total number of output at the last round 
of game is -I- 4. If every sub-game is different from an- 
other, the dimension of probability variables is n^. How- 
ever we have only one equation to solve, G{t) > 0. Even 
if all the sub-games are losing, the solution space is a 
dimensional cubic, each edge of the cubic has a length of 
[0, 1/2). To solve one inequality equation of unknown 
variables is a trivial problem. But if all the subgames 
labeled by 1 are the same, so does the game with label 2. 
We have only two variables which are distributed on the 
whole lattice, it is an interesting question to ask wether 
the combined game wins or loses. 

To find out the final output of flipper network has a 
lot of things in common with many body physics. Espe- 
cially when the flipper are quantum processors, we are es- 
sentially manipulating various entangled quantum states 
from local flipper to reach a global winning output. Many 
coupled identical flippers on a regular lattice can be theo- 
retically summarized into a many body Hamiltonian sys- 
tem. Each flipper represents one losing game, this many 
body system represent a complex game of many coupled 
losing sub-games. We interested in the existence of a 
winning combined game out of many losing games, or 
the opposite case. Although this is no longer a paradox 
due to the strong coupling between neighboring losing 
games, it is still an interesting question. I will take many 
Ising spins on square lattice as an example to show how 
to model a quantum many body system into a combined 
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game of many losing sub-games. We use Pauli matrix 
erf to express a flipper. Each flipper occupied one lattice 
site. The nearest neighboring flippers couple with each 
other to work together. The output of the flipper pair 
afffj for the input state \siSj) is 

h,j = (s.sjla^a'^ls.s^), (50) 

Si = ±1 is Ising spin. The output operator of the com- 
bined game on the whole lattice is Hamiltonian 

H^J2 J^j^^^h (51) 

where < ij > denotes the nearest neighboring lattice 
sites. Jij is the coupling strength. The input capital is 
the difference between the total number of spin -1-1 and 
the total number of spin down —1 in any given initial 
state 

lAnput) = \S1S2S3 ■ ■ ■ Sn) ■ (52) 

For instance, \ipinput) — |1- I7 ■ • • 1 I,)- The 

game rule is to minimize the total energy to reach a 
ground state. Playing one round of game is performed 
by operating the Hamiltonian on the initial configuration 
once. The output state after k rounds of game is 

IV'Lput) =^1^mp«t). (53) 

The spin configuration at each round of game must make 
sure that the total energy E = {tp\H\ip) decreases com- 
paring with last round of game. We take st = +1 as 
the positive capital for the game, if it becomes —1 in the 
end, the game lost. The total game covers many sub- 
games. If we input less {si — +1}, but end up with 
more {si = +1} than our input, the games wins. The 
Hamiltonian Eq. (I51|) is equivalent to the familiar Ising 
Hamiltonian H — J2<ij> '^ij'^i'^j- Theoretical physicist 
always search for the stable spin configurations of ground 
state starting from an arbitrary initial configuration, es- 
pecially in quantum monte Carlo simulations. The so- 
called game paradox appears almost everywhere in many 
different quantum models. In mind of the great inter- 
est on time series in stochastic dynamics, what quantum 
Monte Carlo simulation has abandoned during the pro- 
cess of finding ground states maybe can help people to de- 
sign winning games by combining a large number of lost 
sub-games. On the other hand, the so-called game para- 
dox remind us that designing different non-equilibrium 
game process maybe leads to different physical proper- 
ties of many body system. 

VII. SUMMARY 

When we divide a complex winning game into many 
sections, each section has an input and output, we call 
it a sub-game. Every sub-game must be tested to see 
wether it is a losing or winning game. It is not against 



intuition that if all sub-games are lost. More over, if we 
combine these sub-games together but in a different or- 
der as the original one, it is not against intuition either 
if the combined game is a losing game. No mater it is di- 
viding or combining, we are dealing with the correlations 
between those sub-games. A real game paradox only ex- 
ist when there is no any correlation between sub-games. 
In that case, any conclusion about the total game based 
on the result of sub-games are not reasonable. 

Both Parrondo-Harmer- Abbott game and Parrondo's 
game introduced the correlation of probabilities between 
the nearest two rounds of game. In mind of a fact that the 
first round and the second round are the nearest neigh- 
bor, the second round and the third round are also near- 
est neighbor, and so forth, it forms a long chain of many 
coupled sub-games in history. In my optical analogy, 
the photon propagates through a long path connecting 
many beam splitters. The photon is either reflected or 
transmitted at each beam splitter. The final intensity 
of this photon is the product of the reflection or trans- 
mission coefficient of every beam splitter along this path. 
The more rounds of game one played, the more paths one 
would have. We sum up the intensity of all possible paths 
to reach the final output. If we have more transmitted 
photon than reflected photons, the final gain would be 
positive, the game wins, otherwise the game lost. If the 
photon is reflected for odd number of times, its contri- 
bution to the final gain is negative. If the photon is re- 
flected even number of times, its contribution is positive. 
Since reflection photon can transmit through or reflected 
by other beam splitter, the lost of one local sub-game 
can become positive when a negative photon is reflected 
again by another beam splitter. Both Parrondo-Harmer- 
Abbott game and Parrondo's game can be mapped ex- 
actly into the propagation of a photon through the array 
of beam splitters. 

To implement Parrondo-Harmer- Abbott game, we can 
use the same type of physical beam splitters with dif- 
ferent transmission and reflection coefficient. While to 
implement Parrondo's old game, we need two different 
types of beam splitter: one add a phase shift to its 
reflection, the other does not. There is modular oper- 
ation in Parrondo's old game that the capital in game 
must be checked to see if it can be divided by a integer 
or not. This modular operation is equivalent to a special 
beam splitter from the point view of probability distri- 
bution. If we put this special beam splitter aside and 
only compare the optical diagram structure of Parrondo- 
Harmer- Abbott game and Parrondo's old game, they are 
essentially same. Parrondo-Harmer-Abbott game and 
Parrondo's game simulate the random sequence of comb- 
ing game A and game B, the final game is winning. In 
fact, what really matters is the coupling between the two 
steps in game B, no matter how random the sequence of 
ABAABBA ■ ■ ■ are, it does not break the internal cou- 
pling inside game B. It is through the coupled two steps, 
or in other words, through two neighboring beam split- 
ters, the probability of negative capital is transformed 
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into the probabihty of positive capitaL This hes in the 
heart of a winning combined game out of losing games. 

We come to the conclusion: there is no paradox in 
both Parrondo-Harmer- Abbott game and Parrondo's old 
game, these two games are both reasonable and share the 
same strategy of coupling the nearest two rounds of game. 
This is a good news to gamblers for they can win from 
loss by reasonable calculations without worrying about 
paradox. 

One can design much more complex games by draw- 
ing optical diagrams following the same strategy I used. 
A three-step history-dependent game is designed as an 
example to show how to combine many games. In this 
game, even if all the sub-games lost, the output of this 
combined game can win, lose or oscillating between loss 
and win. Different probability distribution leads to dif- 
ferent results. To calculate the final output, one can draw 
all the paths connecting beam splitters. The final gain 
of each path is derived by multiplying the value of prob- 
ability of all the bonds along this path. Summing up all 
the paths gives the final gain. A real photon propagating 
across array of beam splitters provide a practical way to 
test the complex design. 

Since correlation can make losing games win, it is a 
good strategy to implant more correlations among sub- 
games. I designed a history entangled game with only 
two sub-games. The input state of each sub-game de- 
pends on the symmetric or antisymmetric combination 
of the historical states two steps earlier. In the symmet- 
ric states, two winning sub-games may leads to a lost 
combined game. If one of them wins, the game can win. 
In the antisymmetric state, two losing game results in a 
winning combined game. If one wins, the other lost, the 
output is either lose or win, it depends on specific param- 
eters. In some parameter region, the final gain oscillates 
between win and loss. 

We can implant strong correlations into the capital 
instead of sub-games. This is in case if one can not 
modify the game rule, but one can change the way of 
playing the capital in hand. Suppose one has 1000 
points in the beginning, we cut it into many pieces, 
100 -I- 350 -MO -I- 50 H = 1000. Then we invest the 



capital in game according to a wave pattern, such as 
cos[t]^ -I- sin[t + t]'^, t is the delay time. When this capi- 
tal wave propagates through the lattice of many coupled 
sub-games, the outputs of sub-games are strongly entan- 
gled. To determine the final output of these quantum 
entangled states is equivalent to solve a quantum many 
body physics problem. For example, when we search for 
the spin configuration of ground state on a lattice, usu- 
ally it starts from a random initial spin configuration, 
this initial configuration can be viewed as input capital 
of a game. Then we update the initial configuration by 
applying a functional of Hamiltonian operator, this is a 
process of combined game, every local operator on single 
lattice is a sub-game. The final output of the game is sta- 
ble spin configuration. This is how quantum Monte Carlo 
simulation runs everyday. In fact, Monte Carlo method, 
as its name suggested, comes from gamble games. 

Coupled many sub-games is useful for designing var- 
ious complex Brownian ratchet. The one dimensional 
periodic asymmetric potential is equivalent to a chain of 
many beam splitters. If a flashing Brownian can cover 
two or more potential well within one period, we must 
take into account of the coupling between two or more 
neighboring potential wells. The forward velocity is de- 
termined by summing up continuous two steps of pass 
or refelection. The backward velocity is determined by 
the joint probability of one reflection followed one pass 
or vice verse. The optical model provide a tool to an- 
alyze the relationship between different individual po- 
tential wells. Especially when there are many different 
individual local potential wells distributed along the long 
chain. If each local potential well has different reflection 
and transmission, it is hard to see the physics directly 
from Fokkcr-Planck equation or Smoluchowski equation. 
In that case, we take every local potential well as a dif- 
ferent beam splitter, draw the optical diagram, then one 
will get at least a phenomenological understanding. 
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